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1 Introduction 



In relativistic nuclear (A-A) collisions we search for a macroscopic, deconfined space time 
region in which soft quark and gluon degrees of freedom are present, the so-called quark-gluon 
plasma (QGP) — 'plasma' since it involves freely moving color charges within the limited 
reaction volume. Existence of this phase of matter seems to be an inescapable consequence of 
our theoretical knowledge about the fundamental hadronic interactions, qualitatively rooted 
in quantum chromo-dynamics (QCD). Nevertheless, it is important that we can demonstrate 
QGP formation experimentally, confirming the theoretical paradigm and permitting the 
experimental study of the strong interactions vacuum state. 

Because of its rather limited (nuclear size) volume and hence very short lifespan of 
about 4 fm/c, the study of the QGP phase is a formidable task. Moreover it is far from clear 
that this new phase can be formed and studied with the presently available nuclear beams, 
up to 200A GeV energy at CERN-SPS. However, we have shown in our recent work |Q ^ 
that strange particle production results point towards formation of the deconfined phase. In 
this paper we develop in full the method that involves as the observable strange antibaryons 

^ . Our primary interest is oriented towards (multiply) strange antibaryons for a number 
of reasons: the conventional backgrounds are small since the multi-step processes required 
occur relatively rarely in p-p interactions. Consequently, the relatively high production rate 
of multiply strange antibaryons in A-A reactions, and the central (in rapidity) spectral 
distribution, are indicative of a 'collective' formation mechanism. In the QGP reaction 
picture it is the ready made high density of (anti) strange quarks which leads to highly 
anomalous yields of multiply strange particles. 

The use of specific hadronic particles in the search for QGP is only possible if the evolution 
of the final state is such that information about key properties of the primordial source is 
retained during the formation and evolution of the final state hadrons. This is observed 
experimentally indicating that there is no opportunity for the strange antibaryons emerging 
from deconfined state to re-equilibrate (annihilate) 0, thus implying a sudden transition 
into the final hadronic particles Of course it is far from certain today that QGP is 

formed already in the nuclear collisions below ^/s = lOA GeV. Several studies addressed 
recently the question [0, |^, ||, [10| how the current data can distinguish between the 
formation of a central fireball consisting of a dense and very hot hadron (resonance) gas and 
the creation of a transient deconfined state such as the quark-gluon plasma (QGP) which 
subsequently hadronises. We note that even if this final state hadronisation involves a full 
re-equilibration, then the memory of the QGP should in some circumstance not be fully 
lost, since the deconfined state is often richer in entropy than the confined state 0, |Tl|, |2^, 



13], which becomes visible as particle excess in the final state. Today, based on diverse 



results obtained at 200A GeV [g] it appears that we have a very favorable situation with the 
observable particles very much representing the expected properties of the primordial QGP 
state, while the results from BNL-AGS at 10-15A GeV are indicating that if QGP was formed. 



it has undergone a full re-equilibration transition []T4[. Consequently, somewhere between 
the two energy regimes a major change in the reaction mechanism occurs. It appears to be 
of considerable importance to find in which energy range this change occurs and how other 
observables vary with energy in this changeover region. 
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Since the QGP picture works so well at 200A GeV, it is natural to expect that its 
formation is also occurring at somewhat lower energies. Therefore, within a dynamical 
collision model, using the equations of state (EoS) of QGP, we study in this paper how the 
thermal and chemical properties of the fireball of the dense hadronic matter change with 
energy. We predict here the production of strange particles as function of energy. Our 
primary goal is to find how the strange particle abundance in the final state is related to the 
conditions prevailing in the early moments of nuclear collisions. We develop an approach 
in which the collision force is counteracted by the thermal radiative pressure. The tacit 
assumption that we make is that the thermal equilibrium of quarks and gluons arises rapidly 
as a result of the primary interactions. This introduces without an explicit mechanism almost 
70% of the final state entropy []r^ — this is the excess QGP entropy visible as the excess 
multiplicity of produced particles. Thus in our approach we do not resolve the important 
issue, what is the energy and more generally the conditions required in the collision for the 
formation of the deconfined state. On the other hand, the principal virtue of the thermal 
model framework adopted here is that the spectra and particle abundances can be described 
in terms of a few parameters which have simple physical meaning. 

Aside of the rather demanding multi-strange antibaryon signature of the QGP, there is 
the more readily available, but less specific overall abundance of strangeness. Strange particle 
signatures of dense baryon rich hadronic matter have been theoretically considered for nearly 



two decades [|T6|. Because of the intrinsic difficulties related to the detection of strange 
particles within the large multiplicity of hadrons produced in relativistic nuclear collisions 
only recently results have been presented which allow a test of the different theoretical models 
(for a recent review of available results, see Ref. fl^l)- Because strangeness production 
processes constitute a bottleneck in chemical equilibration of conventional confined strongly 



interacting matter [18, M, strangeness production in QGP has been widely studied 



1^, |T9|, ^ . In the theoretical description several microscopic features combine to produce this 
enhancement: there is a high quark-gluon density and the production threshold of strange 
quark pairs in QGP is below the mean energy (temperature) of the constituents present. 
The computed strangeness relaxation time constant Tg < 6 fm/c towards the attainment 
of the strange phase space equilibrium abundance in the QGP phase is of the magnitude 
similar to the estimated lifespan tqgp — R/c of the dense QGP state, where R is the 
characteristic size parameter of the QGP phase. Consequently, even if the magnitude of both 
these time constants is still rather uncertain, their similarity is of considerable importance 
for the exploration of the QGP phase: varying the size of the colliding nuclei, or the mean 
impact parameter in collision, we can change the size R of the interaction region and thus 
change the lifespan of the dense state. Since Xg depends primarily on the temperature (mean 
glue energy), we could explore how this impacts the strangeness chemical equilibration and 
compare with the theoretical predictions. 

A crucial confirmation of the reaction picture would arise if the changes observed in the 
chemical equilibration of strangeness would be equally attainable when the energy content 
in collision is varied. More importantly, this could be done at the maximum volume avail- 
able, thus assuring that unwanted dynamical variations such as changes in longitudinal flow 
(transparency) with changes in impact parameter were minimal, while the highest accessible 
particle and energy density is explored in each case. For the strange antibaryon signature 



4 



J. Letessier, J. Rafelski and A. Tounsi 



this approach is much preferred, as the abundance of particles produced, and their centrahty 
remains assured. Since the (enhanced) production of (multi) strange antibaryons is widely 
considered to be more specifically related to the deconfinement than the strangeness enhance- 
ment in general, it would be of considerable importance to obtain data on (multi) strange 
antibaryon production as function of energy. 

Within our approach one of key tacit assumptions is that there is formation at central 
rapidity of an energy and baryon-rich, fireball. We believe that the pattern of diverse particle 
production processes confirms the view that even the 200A GeV collisions are leading to the 
formation of a domain of space-time in which not only energy, but also considerable baryon 
number is present. For example, the rapidity distribution of A and A as measured by the 
NA35 experiment |^ affirms such a picture of the A-A reactions. These results show that 
the S-Ag/W/Pb and even the S-S collisions at 200A GeV are far from the limit of baryon- 
transparency in which the valence quarks of projectile and target are presumed to leave 



the central rapidity region |23]. Relatively large stopping behavior at 200A GeV could be 



expected: assuming that a valence quark carries, in a typical case, about 1/6 of the energy 
of a nucleon, and placing the observer in the CM frame, we find that we would have two 
1.5 GeV partons colliding at CERN-SPS energies — at BNL-AGS energies (10-15A GeV 
) a head-on collision interpreted in this picture would occur between partons of 0.4 GeV. 
Moreover, we have in each hadron a rather wide distribution of parton energies, and thus 
many collisions between sea and valence partons are occurring at still considerably lower 
energies. Just a few elastic or inelastic interactions will suffice to stop many partons in the 
central fireball. It is thus not surprising that in the particle spectra there is little room, if 
any, for the separation in rapidity of energy and baryon number. 

In this consideration it is important to remember that the more easily analyzed and un- 
derstood symmetric A-A collisions at small impact parameters the surface (corona) nucleons 
will undergo interactions which resembles the normal p-p collision environment. These p-p 
like components should not be confused with the more specific elements of the A-A inter- 
actions — strange particles are a good tag on these specific interactions, since strangeness 
production in A-A collisions is enhanced compared to p-p and p-A based expectations, 
and because there are a number of particles (such as A) which are rather rarely produced 
in p-p collisions, in contrast to A-A interactions, and which show, very strongly in these 
interactions, features we would interpret as originating in the central fireball. 

Our paper is organized as follows: in next section ^ we explore the conditions reached in 
the early deconfined QGP fireball formed e.g. in the CERN-SPS collisions. We then compare 
in section |^ the results we obtain with those derived for the S-W/Pb at CERN-SPS and 
Au-Au interactions at BNL-AGS. We obtain an excellent agreement with all the available 
data. We also make predictions about strange particle production for the forthcoming round 
of experiments and about the expected particle yield behavior as the energy of the collision 
is varied. 



2 Properties of the QGP State 
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2.1 The Thermal Approach 

Application of the thermal model to hadronic processes implies that on the prevailing time 
scale of the interaction, thermalization of the relevant degrees of freedom is rapid. This 



means that thermal features are seen in the particle production processes which are not 
rooted directly in the microscopic features of strong interactions (QCD). We can expect that 
the leading particle and flow effects diminish in their significance as the size of the interaction 
region increase. Thus soft hadronic physics in the A-A case should be easier to understand. 
The original motivation to study these large nuclei in highly relativistic collisions was just 
the hope that we will unravel the soft hadronic physics which has remained in veil when 
studied with the p-p interactions. 

It is convenient to incorporate in the model development the particle spectra as function 
of the rapidity y and transverse mass m± of a particle 

E + p, 



2/ = ^ In 1^-^— ^ j , E = m_LCOsh.y, mj_ = ^m^ + p'j_, (1) 

where '_L' is perpendicular to the collision axis While m± is invariant under Lorenz 
transformations along the collision axis, the particle rapidity y is additive, that is it changes 
by the constant value of the transformation for all particles. This allows to choose the 
suitable (CM — center of momentum) reference frame characterized by its rapidity ycu for 
the study of the particle spectra. 

The 'thermal' hypothesis presupposes that the shape of these spectra is sufficiently similar 
to allow a reduction of all data to just one 'Boltzmann' spectral shape centered around the 
fireball rapidity ycu- 

'£L = iV,e-^'''/^ = NiC- --My-ycM)m*JT ^2) 
d?p 



The parameters of each particle distribution [25| include the inverse slope T ('temperature 



of the distribution, centered around the ycu- In particular the temperature T derived 
from the transverse mass distribution should be common for all types of particles produced 
by the same thermal mechanism. This relates many particle spectra to each other and 
reduces the number of observables considerably. Similarly, |/cm should be also the same for 
different particles originating in the same thermal source. In addition, for each particle there 
is the normalization constant Ni. We will discuss below how these constants are determined 
by a chemical equilibrium parameters, aside of the volume V . Let us note in passing that 
the presence of a common inverse slope parameter for different particle spectra, which is 
different from the one found in p-p collisions is well established experimentally. For example, 
for the 200 A GeV S-W/Pb interactions the inverse slope parameters of strange baryons and 
antibaryons and for example of high m_L neutral mesons |^ are all consistent with 
T = 232 MeV. 

Since a wealth of experimental data can be described with just a few model param- 
eters, this leaves within the thermal model a considerable predictive power and a strong 
check of the internal consistency of thermal approach. Specifically, in the directly hadro- 
nising off-equilibrium QGP-fireball model there are 4 parameters (aside of T and ycm) 
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characterizing all particle spectra: two fugacities Aq, As, and two particle abundance non- 
equilibrium parameters, the strangeness occupancy factor 7s and the ratio Rl related to 
meson and baryon abundances (see section 3^). Moreover, these physical parameters can 
be determined using a dynamical picture of the collision, in which the input is derived from 
more general qualitative conditions of the colliding system, such as the energy content or 
stopping power. Thus the validity of thermal and chemical equilibrium can be conclusively 
tested, comparing the predicted particle spectra and yields with theoretical predictions with- 
out the need or capability to modify and adapt the description to each new result which is 
reported. 



2.2 Parameters of the Fireball 

We now introduce and discuss in qualitative terms the model parameters and their temporal 
evolution during nuclear collision. Practically all results we will obtain will hinge on the 
time evolution scenario we adopt. At the present time the experimental constraints are not 
limiting our reaction and evolution picture in a unique way. The picture we adopt is in 
qualitative, and as we shall see, also in quantitative agreement with a number of important 
experimental results, and it is in accord with the general properties of the strong interactions 
and hadronic structure widely known and accepted today. 

2.2.1 Temporal Evolution of the Fireball 

We suppose that the relevant time development stages of the relativistic nuclear collision 
comprise: 

1. The pre-thermal stage lasting perhaps 0.2-0.4 fm/c, during which the thermalization 
of the initial quark-gluon distributions occur. During this time most of the entropy 
obtained in the collision must be created by mechanisms that are not yet understood 
— this is also alluded to as the period of de-coherence of the quantum collision system. 
Our lack of understanding of this stage will not impact our results, as the reason that 
we lack in understanding is that the hadronic interactions erase the memory of this 
primordial stage. 

2. The subsequent inter-penetration of the projectile and the target lasting about ~ 1 
fm/c, probably also corresponding to the time required to reach chemical equilibrium 
of gluons G and light non-strange quarks q = u, d . 

3. A third time period (3-5 fm/c) during which the production and chemical equilibration 
of strange quarks takes place. During this stage many of the physical observables 
studied here will be initiated. 

4- Hadronisation of the deconfined state ensues: it is believed that the fireball first ex- 
pands at constant specific entropy per baryon, and that during this evolution or at its 
end it decomposes into the final state hadrons, under certain conditions in a (explosive) 
process that does not allow for re-equilibration of the final state particles. 
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In the sudden hadronisation picture of the QGP fireball suggested by certain features seen in 
the analysis of the strange antibaryon abundances for the 200A GeV nuclear collision data 0, 
^ |], the hadronic observables we study are not overly sensitive to the details of stage 4- Akin 
to the processes of direct emission, in which strange particles are made in recombination- 
fragmentation processes 0, the chemical conditions prevailing in the deconfined phase are 
determining many relative final particle yields. If the hadronisation occurs as suggested by 
recent lattice results at a relatively low temperature (e.g. 150 MeV), the total meson 
abundance which is governed by the entropy content at freeze-out of the particles, is found 
about 100% above the hadronic gas equilibrium expectations |Tl|. This is consistent with 



the source of these particles being the QGP ITT], The freeze-out entropy originates at 
early time in collision since aside of strangeness production which is responsible for about 
10% additional entropy there is no significant entropy production after the initial state has 
occurred 0]. Similarly, the relatively small thermal abundance of baryons must be enhanced 
by the factor 5 in order to maintain a ratio Rl = 0.4 [Q]. 

On the other hand the experimental results obtained at 15A GeV are consistent with the 
thermal equilibrium hadronic gas state expectations [|14|. In the event that the source of 



these particles should indeed be a QGP fireball, a slow re-equilibration transition should be 
envisaged here. In this scenario the details of hadronisation mechanisms are of even lesser 
relevance since the equilibrium state is reached among the final hadron gas particles. 



2.2.2 Temperature 

During the initial contact of the two nuclei, as soon as stage 1 is reached, it is likely that 
the most extreme conditions (highest temperature) prevails. Subsequently, temperature de- 
creases rapidly, primarily due to two mechanisms: there is rapid quark and gluon production, 
leading towards the chemical equilibration, and second, the fireball expansion/emission cool- 
ing occurs. We will account for these two effects balancing energy per baryon or entropy 
per baryon of the fireball as appropriate for the evolution stage. In this way, we implement 
during the stage 2 a one dimensional, hydrodynamical evolution of the dense deconfined 
hadronic matter, in addition allowing for the approach to chemical equilibrium. We en- 
counter a considerable drop in temperature between the initial stage 1 and the final state 
4 freeze-out point, but the entropy content which determines the final particle multiplici- 
ties remains nearly constant: aside of the initial state entropy formation all the additional 
increase is due to the formation of the strangeness flavor. 

The different stages of the evolution are characterized by the temperatures: 

Tth temperature associated with initial thermal equilibrium, 

I production of q, q, G, expansion; 

Tch temperature of chemical equilibrium for non-strange quarks and gluons, 

I production of s, s quarks and fireball expansion; 

Tq temperature of maximal chemical equilibrium: 'visible' temperature, 

I fireball expansion/particle radiation; 

Tf s temperature at freeze-out for non-strange and strange particles, 

with obviously Tth > Tch > Tq > Tf. 
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In the transverse mass spectra of strange (anti)baryons a temperature T± is found. If 
the final state particles emerge directly without re-equilibration from the fireball ^ , this 
observed temperature (Tj_ = 232 ± 5 in S-A collisions at 200A GeV) in the particle spectra 
would be closely related to the full chemical equilibration temperature Tq : subsequent to 
the establishment of the conditions at To we have either directly the emission of particles 
and thus we have T± < Tq, or there is collective, (so called transverse) radial flow in the 
hot matter, in which fraction of the thermal energy is converted into the flow energy. When 
the final state particles emerge from the flowing surface, they are blue-shifted by the flow 
velocity. This Doppler shift effect restores the high apparent Tq in high m± particle spectra 

Tx<To^./^+^Tf, T^<To = ./l±^T,; (3) 

— strange and non-strange particles may not originate in exactly the same condition and 
hence we introduced here the option to have different freeze-out temperature. In view of the 
smaller strange particle cross sections we expect Tg > Tf. 

As suggested in Eq. (H), the temperature Tq at the onset of flow is close to the inverse 
slope spectral temperature implied by the high ■m± (strange antibaryon) particle spectra. 
Despite our ignorance of the freeze-out mechanisms and conditions, we believe that the 
uncertainty in the value of the initial temperature as derived from the value of T± in Eq.|^ 
is not large. If QGP phase is directly dissociating by particle emission, this is trivially so. 
If there were to be substantial flow, one can assume some temperature Tq, and given the 
EoS, compute accurately the hydrodynamic radial expansion at the high ■m± ~ 2 GeV 
considered here, using Eq. (^ we find that the inverse slope T± of the particle distributions 
is equal or a bit smaller than Tq. 



2.2.3 Fireball Rapidity 

The fireball is created in central symmetric collisions at the CM-rapidity of the N-N system, 
which is for relativistic systems just is 1/2 of the projectile rapidity. For asymmetric collisions 
such as S-Au/W/Pb the CM rapidity depends on the ratio of the participating masses of 
the projectile Ap and target At nuclei. The center of momentum frame is |^ (neglecting 
small corrections): 

l/CM = y--ln-^. (4) 

Assuming small impact parameter collisions with a suitable central trigger, all projectile 
nucleons participate while the target participants A^ can be estimated from a geometric 



'interacting tube' model [^. This model reproduces well the central value of rapidity 
center particle spectra in the specific case of 200A GeV S-Au/W/Pb interactions. One in 
particular finds ycu = 2.6±0.1, the uncertainty arising from the impact parameter averaging 
and variations of the surface nucleon participation. Once the central rapidity is defined, and 
the ratio of participating projectile and target masses is known, it is possible to determine 
the CM-energy involved in the interaction. It is from this simple kinematic considerations 
that we derive the energy values presented below in section p.4|. 
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2.2.4 Particle Fugacities 

Among the chemical (particle abundance) parameters there are the well known particle 
fugacities, which allow to conserve flavor quantum numbers. Three fugacities are introduced 
since the flavors u, d, and s are separately conserved on the time scale of hadronic collisions 
and can only be produced or annihilated in particle-antiparticle pair production processes]]. 

The fugacity of each hadronic particle species is the product of the valence quark fugac- 
ities, thus, for example, the hyperons have the fugacity Ay = Xu^dK- Fugacities are related 
to the chemical potentials /ij by: 

X - = e^'/'^ , = A^^ i = u, d, s. (5) 

Therefore, the chemical potentials for particles and antiparticles are opposite to each other, 
provided that there is complete chemical equilibrium, and if not, that the deviation from 
the full phase space occupancy is accounted for by introducing a non-equilibrium chemical 
parameter 7 (see below). 



2.2.5 Phase Space Occupancy 

Thermal and chemical equilibria are two very different phenomena. In general, the produc- 
tion of particles is a considerably slower process than elastic hadronic collisions, and thus 
even if we assume in our work a thermal equilibrium scheme, we should not expect the 
chemical equilibrium to be present. In addition, there is the relative and absolute chemical 
equilibrium. In the former, the particle abundances are in relative equilibrium with each 
other, in the latter the total particle yields are just fllling the full available phase space — 
relative chemical equilibrium is in general easier to attain than the absolute chemical equi- 
librium. Our picture of rapid QGP flreball disintegration in which no equilibration takes 
place in the flnal state, implies that chemical (abundance) non-equilibrium features should 
be present in the flnal state. 



Strangeness flavor provides an interesting example of the above. Calculations |T^, |2T| of 
the chemical relaxation constant show that in general it will not fully saturated the available 
phase-space. Therefore, we need to introduce the associated off-equilibrium parameter 7s. 
Since the thermal equilibrium is, as discussed above, established within a considerably shorter 
time scale than the (absolute) chemical equilibration of strangeness, we can characterize the 
saturation of the strangeness phase space by an average over the momentum distribution: 

^ J d^d^xns{p,x;t) 
~ Jd^db:n^{p,x) ' 

where is the equilibrium particle density. The factor 7s thus enters the phase space 
momentum Boltzmann distribution as a multiplicative factor, and with the x dependence 
contained solely in the statistical parameters we have: 

nM X- 1) ~ i,{t)n'^{p- T{x, t),fi,{x, t)) . (7) 



"'^In many applications it is sufficient to combine the light quarks into one fugacity Aq = AdAu , fiq = 
{flu + Md)/2. The slight asymmetry in the number of u and d quarks is described by the small quantity 
(5/i = /id ^ 7 which may be estimated by theoretical considerations 
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A further refinement sliould be noted: wlien tlie quantum aspects of the particle distribu- 
tions are incorporated and the maximum entropy state of an isolated physical system (closed 
system) is obtained, the absolute chemical equilibrium coefficients 7^ enter as multiplicative 
coefficients in front of the Boltzmann factor within the quantum Bose/Fermi distribution 



15 



B,F 1 /o\ 

= T 1 ■ 

along with the fugacity factors . 

The final freeze-out value of 7s is a parameter which impacts the particle ratios in a 
two fold way: there is the simple proportionality (in the Boltzmann limit) of the particle 
yields to 7", — where n is the number of strange quarks contained in the hadron consid- 
ered, and furthermore, the value of 7s impacts the cooling that the fireball undergoes to 
produce strangeness. Unexpectedly, it is rather straightforward to extract from the strange 
antibaryon experimental particle yields [^, |[ |^, [10| the value of 7s. Roughly speaking, 7s 



makes its appearance in all particle ratios in which we compare the abundances involving 
different strangeness content. The measurement of 7s is thus an important step towards the 
understanding of the behavior of highly excited hadronic phases: 7s can be studied varying 
a number of parameters of the collision, such as the volume occupied by the fireball (varying 
size of the colliding nuclei and impact parameter), the trigger condition (e.g. the inelas- 
ticity), the energy of colliding nuclei when searching for the threshold energy of abundant 
strangeness formation. 

The theoretical dynamical model to investigate 7s (i) has been developed to considerable 



detail [32, 33] — it arises from a standard population evolution equation and its interplay 
with the expansion-dilution-dissociation of the fireball. For a QGP state with its fast gluonic 
strangeness production even the natural short fireball lifetime of only a few fm/c should be 
nearly sufficient to reach values of 7s — 1. Detailed balance assures that the production and 
annihilation processes are balancing each other as 7s — > 1 . At large times the approach to 
equilibrium takes the form 

l-7s^e-*/"% t>rs, (9) 
where Tg is the relaxation time for strangeness equilibration which can be computed using 



standard QCD methods. A recent evaluation of this work |21| confirmed that the glue fusion 



processes |TB| are dominating the strangeness production rates in QGP, with Tg — 2 fm/c 
for the here relevant T ~ 250 MeV temperature range. The experimental results yield 
7s = 0.7-1, imply according to Eq. (^ a lifespan of the the plasma state to be greater than 
3 fm. 

In a purely HQ fireball, with its expected much longer strangeness saturation time scale, 
small values of 7s ~ 0.1 like those extracted from A^-A^ collisions should prevail [^. A 
measurement of 7s thus in principle provides important information on the strangeness pro- 
duction time scale and hence a relatively large value of 7s requires some new physics feature 
in the structure of the collision fireball, which we like to associate with QGP formation. 
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2.3 QGP Equations of State 

The partition function of the interacting quark-gluon phase can be written as: 

9iias)V 

iGQGP 



lnZQGP= J2 ^^^^/±ln(l±7iV-v/^^^?(^ (10) 



where i = G, q, q, s, s, with Aj = Aj ^ and = ji. We take into account the QCD 
interaction between quarks and gluons by allowing for thermal masses 

m^{T) = {m'^r + {cTf . (11) 

For the current quark masses we take: 

= 5 MeV, m° = 160 MeV, = . 

We have cx as, as being the QCD coupling constant. Considering the theoretical uncer- 
tainty regarding the coefficient c(as), we explored two different approaches: 

1. We fix c = 2, arising for a^ ~ 1 (the exact value was not of essence), and allowing 
for another effect of the QCD-interaction, the reduction of the number of effectively 
available degrees of freedom. We implement the following effective counting of gluon 
and quark degrees of freedom, motivated by the perturbative QCD formula: 

/ 15a< 
gc =16 5(G(as) = 16 1 



4:71 



( 50a, , 

^7,_T = 6 ^ (?i_T(a,) = 6 (l - ^ ) , (12) 

( 

gi^B =6 gi-B{as) =61- 2 — 

V TT 

where i = u, d, s. In Eq. (0) two factors are needed for quarks: the factor gi^T 
controls the expression when all chemical potentials vanish (the term in the partition 
function for massless quarks) while gi-B is taken as coefficient of the additional terms 
which arise in presence of chemical potentials. We took a^ = 0.6 which turned out to 
be the value best suited for the experimental data pointsQ 

2. A fit proposed by Kampfer et al. 0, ^ made to the lattice data for the SU(3) gauge 
sector and extended to quarks 

= 3.8a^(T); (13) 

471 

' ^ ' 9 ln(T/T, + 0.023)2 ' ^ ^ 

where the form of a™ is motivated by aslq"^). The critical temperature of the QGP 
transition to confined phase is chosen at = 150 MeV and the infrared cutoff at 
0.023T,. 



^This value arises at T = 2.2 Tc in the fit p% pq| to lattice gauge data without quarks. 
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Figure 1: Properties of different QGP-EoS: constraint between temperature T and light quark 
fugacity Aq for a given energy content per baryon E/B = 4.3 (thin hues) and 9.6 GeV (thick hues). 
Upper sohd line: procedure 1 of EoS (see text), lower solid lines joining onto the long-dashed lines 
(second solution branch): procedure 2 EoS and short-dotted: free quark-gluon gas. 



The main difference between these two approaches is how the quark and gluon degrees of 
freedom are suppressed near to the critical temperature. In both cases these phenomeno- 
logically motivated procedures are a bold extension of the established perturbative and/or 
lattice QCD results. 

Given these interaction effects, the Bose/Fermi distributions for quarks/gluons, Eq. (p!0|), 
allow us to obtain in a consistent way any physical property of this model EoS of the QGP 
phase. These are shown in Fig. |^. We show here in T-Xq domain lines of fixed energy 
per baryon (thick lines: 9.6 GeV, thin lines 4.3 GeV) — solid lines are for the procedure 
1, dashed lines are the ideal gas (no interaction) lines, long dashed is the second solution 
branch of the procedure 2, which attaches to the primary solution branch also shown as 
solid lines. Note that the second solution branch is probably an artifact, as it tangents the 
critical temperature, where this method is questionable. Interestingly, we note that at high 
T the perturbative result is just below the results of procedure 2, indicating that the change 
in energy density and baryon number density due to thermal masses is equal, and thus the 
ratio is similar to the free gas result. There is substantial difference between the procedure 
1 and the free gas, which originates in different changes of the degrees of freedom due to 
perturbative interactions in the dominant terms of energy density and baryon density. 
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This representation of the ratio of two densities amphfies the impact of interactions 
and shows how these are incorporated in the different approaches. We see that while in 
the procedure 2 fitted to lattice data the treatment of the energy density is particularly 
accurate for baryon-less plasma, for finite baryon density of interest here the results are 
of definitively questionable validity. In particular, near to the critical temperature (taken 
here to be = 150 MeV) the large thermal mass of quarks eliminates the baryon density 
effectively and thus it is not even possible to find for a given ratio of energy density to baryon 
density a reasonable solution for T ~ Tc while for 'high' T essentially the free quark-gluon 
gas result is reproduced. 

Since our present work is primarily concerned with the baryon density dependence of 
the diverse observables near to Tc, presumably < T < 2Tc and /ig ~ T, it is clear that 
we cannot use the model 2 without prior major theoretical improvements applicable to high 
baryon density domain, such as is the inclusion of the chemical potential dependence in 
and in the thermal masses, and for fermions, the effect of Fermi blocking on the polarization 
function. We thus assume the procedure 1 with parameters as described above. 

2.4 Energy Content in the Fireball 

We now commence the evaluation of the initial conditions reached in the collision: an impor- 
tant constraint arises from the energy per baryon content in the fireball. The considerations 
of the previous section allow to relate a given energy per baryon]^ E/B to the statistical 
parameters T and Aq. On the other hand, the collision energy gives 

E tjeEcm Ecu 

^ 

B ^^B^part ^part 

where Apart is the number of nucleons participating in the reaction. The last equality follows 
when the stopping fractions are equal — the experimental particle spectra we are addressing 
here, and in particular the visible presence of baryons in the central rapidity region, are 
implying that this is a reasonable assumption for the current experimental domain. In 
consequence, the energy per baryon in the fireball is to be taken as being equal to the 
kinematic energy available in the collision. In the current experiments we have the following 
kinematic energy content: 



Au-Au at 10. 5 A GeV - 


E/B = 


2.3 GeV , 


Si-Au at 14. 6A GeV 


E/B = 


2.6 GeV, 


A-A at 40A GeV 


E/B = 


4.3 GeV, 


Pb-Pb at 158 A GeV 


E/B = 


8.6 GeV, 


S-W/Pb at 200A GeV - 


E/B = 


8.8 GeV , 


S-S at 200A GeV 


E/B = 


9.6 GeV , 



Note that above we assumed collision with the geometric target tube of matter when the 
projectile is smaller than the target, see section p.2| . The specific energy content E/B, given 
EoS, establishes a constraint between the thermal parameters. In Fig. ^ we show in the T-Aq 

■^Here B is the baryon number. To avoid confusion the bag constant is denoted B. 



(15) 
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Figure 2: QGP-EoS according to procedure 1 constraint between temperature T and light quark 
fugacity Aq for a given fireball energy content per baryon EjB appropriate for the BNL-AGS and 
CERN-SPS collision systems. Left to right: 2.3 (Au-Au), 2.6 (Si-Au), 4.3 (A-A), 8.6 (Pb-Pb), 8.8 
(S-PB/W) and 9.6 (S-S) GeV. See text for a discussion of experimental points. 

plane the lines corresponding to the constraint on the QGP-EoS according to procedure i, 
see section ^.3| , arising from fixing the energy per baryon at kinematic value listed above 
(rising from bottom right to left). In the middle we show the lowest CERN-SPS accessible 
energy, 4.3 GeV, which bridges the current CERN-SPS domain shown to the left to the BNL 
region on the lower right. The experimental crosses show the values of Aq arising in our data 
analysis 0, 14], combined with the inverse slopes temperatures, see section 2^, extracted 



from transverse mass particle spectra. The fact that the experimental results fall on the 
lines shown in Fig.^ is primarily due to the choice = 0.6 — as this is the usual value in 
this regime of energy it implies for a QGP fireball EoS hypothesis that the assumption that 
stopping of energy and baryon number is similar deserves further consideration. 



2.5 Pressure Balance 

We have seen that the QGP energy per baryon constraint between Aq and T (see Fig.^) allows 
for a good agreement with experiment There remains the issue what physical constraint 
or principle determines which of the possible pair of T, Aq values along the individual curves 
in Fig.|^ is experimentally recorded by the cross shown. We have explored the properties 
of the QGP phase along these lines of constant energy per baryon and have noticed that 



Strange Particles Yields from QGP 



15 



with increasing T the pressure in the QGP phase increases, and that the experimental points 
coincide with the dynamical pressure generated in the collision. This confirms the intuitive 
idea that the initial conditions reached in the central fireball arise from the equilibrium 
between the fireball internal thermal and external compression pressure. 
This condition takes the form IST 



Pth(r,Ai,7,) =Pdyn + Pvac. (16) 

The thermal pressure follows in usual way from the partition function 

P,^ = T/V\nZ, (17) 

where aside of the temperature T, we encounter the different (quark and gluon) fugacities 
Aj and the chemical saturation factors 7j for each particle. For the vacuum pressure we will 
use: 

P,ac = B- O.lGeV/fm^ . (18) 
The pressure due to kinetic motion follows from well-established principles, and can be 



directly inferred from the pressure tensor 



T^{x) = JpVf{x,p)d'p, t,j = 1,2,3. (19) 
We take for the phase-space distribution of colliding projectile and target nuclei 

fp,T{x,p) = pp^t{x)6%p±pcm) , (20) 

and hence in Eq. ( PD| ) = ±Pq-^/ Eqm- We assume that the nuclear density is uniform 
within the nuclear size, po = 0.16 /fm^. 

To obtain the pressure exerted by the flow of colliding matter, we consider the pressure 
component with j being the direction of vcm- This gives 

Pdyn = V.P0§^ . (21) 
-C/CM 

Here it is understood that the energy Eqm and the momentum pcM are given in the nucleon- 
nucleon CM frame and rjp is the momentum stopping fraction — only this fraction < r^p < 1 
of the incident CM momentum can be used by a particle incident on the central fireball (the 
balance remains in the un-stopped longitudinal motion) in order to exert dynamical pressure. 
For a target transparent to the incoming flow, there would obviously be no pressure exerted. 
In Fig.^ we show, as function of momentum stopping rjp the collision pressure P. For 
?7p ~ 0.5 the value of the pressure P is close to the value found in the QGP phase (see 
table 1). We note that if indeed in Pb-Pb collisions the momentum stopping increases, the 
dynamical pressure can double, causing a considerable increase in the baryon compression 
and the associated slow rise of Aq as the size of the collision system increases. 
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Figure 3: The collision pressure P as function of momentum stopping r/p for different values of 
E/B — 2.3, 2.6, 4.3, 8.6, 8.8 and 9.6 GeV (from bottom to top, solid line is for 8.8 GeV). 

2.6 Final State Statistical Parameters 

We are here primarily interested in determining the properties of the (near) chemical equilib- 
rium state of the QGP-fireball prior to its disintegration. During the time the nuclei interact 
in collision, we assume rapid thermal equilibration of the fraction t] of the kinetic energy. 
The evolution of particle numbers towards chemical equilibrium occurs in this stage under 
the constraint of a given constant (compression-) pressure and energy per baryon. Given 
values of Eqm and ?7p, as well as taking rj^ = rj^ we have prescribed a particular energy per 
baryon E/B in the fireball and we can concurrently solve the pressure equilibrium constraint 
given by Eq. ([TBI) to determine the pressure. Since the pre-equilibrium QGP phase is not 
observable with hadronic probes, much of the detail of the evolution is lost. As we find, it 
suffices to use the energy and pressure constraints to determine the properties of the fireball 
at the time the nuclear collision terminates. We make here the plausible hypothesis^ that 
when the collision has terminated (at about 1.5 fm/c in the CM frame), the u, d quarks 
and gluons have reached their chemical equilibrium, 7q ^ 1, 7g — *■ 1. Furthermore, the 
strange fiavor is still far from equilibrium and we choose 7s ~ 0.15 appropriate for strange 

''The finite baryon density and baryon number conservation in the fireball force onto the system a rather 
large quark density, which is there from the beginning and needs not be produced; gluons are more easily 
produced than quark pairs and thus presumably their number catches up with the quark number by the 
time the collision has terminated. Note that in baryon-free central region environments expected at much 
higher RHIC/LHC energies, the approach to chemical equilibrium can be different. 
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Table 1: Properties and evolution of different collision systems. 



Phase 


<s - s>= 


E/B [GeV] 


space 


2.6 


4.3 


8.8 


8.6 


8.6 


occiinancv 


A = 1 


f] = 1 


f] = 1 


7] = 0.5 


r/ = 0.75 


rj = 1 






Au-Au 


Pb-Pb 


S-Pb 


Pb-Pb 


Pb-Pb 






0.260 


0.361 


0.410 


0.444 


0.471 


^, no 


\ 

Ag 


9.95 


3.76 


1.78 


1.91 


2.00 




Hg/B 


0.20 


0.54 


1.55 


1.36 


1.25 






3.00 


3.13 


5.12 


3.89 


3.77 






0.00 


0.13 


2.12 


0.89 


0.77 




0.02 


0.06 


0.16 


0.14 


0.13 




^th V / ini J 


0.46 


0.76 


0.79 


1.12 


1.46 


^, n nQ 


Pb 


3.34 


3.30 


1.70 


2.44 


3.18 




S/B 


11.8 


18.8 


40.0 


35.8 


33.4 




Trh fGeVl 


0.212 


0.263 


0.280 


0.304 


0.324 


-y = 1 


A 


4.14 


2.36 


1.49 


1.56 


1.61 




n IB 


0.56 


1.08 


2.50 


2.24 


2.08 




n IB 


3.11 


3.51 


5.16 


4.81 


4.62 


^ = 1 


n-lB 


0.11 


0.51 


2.16 


1.81 


1.62 




n-lB 


0.05 


0.11 


0.25 


0.22 


0.21 




P.h [GeV/fm^l 

cn V / 1111 J 


0.46 


0.76 


0.79 


1.12 


1.46 


o/o = 15 




3.35 


3.31 


1.80 


2.45 


3.19 




S/B 


12.3 


19.7 


41.8 


37.4 


34.9 






1 


1 


0.8 


1 


1 


7.= 1 


To [GeV] 


0.184 


0.215 


0.233 


0.239 


0.255 




\ 


4.14 


2.36 


1.49 


1.56 


1.61 


7.= 1 


njB 


0.56 


1.08 


2.50 


2.25 


2.09 


njB 


3.11 


3.51 


5.12 


4.81 


4.60 




ng/B 


0.11 


0.51 


2.12 


1.81 


1.62 


7. = 0.8 


n-s/B 


0.34 


0.68 


1.27 


1.43 


1.33 


or 


Po [GeV/fm3] 


0.30 


0.41 


0.47 


0.54 


0.71 


7. = 1 


Pb 


2.17 


1.80 


1.05 


1.19 


1.56 




S/B 


14.5 


24.0 


49.5 


46.5 


43.4 



quark relaxation time 7 times larger than the light quark one |21]. Because the QGP phase 



is strangeness neutral we have As = 1. The remaining statistical parameters T^i and Aq 
are now fixed by the EoS and are shown with other interesting properties of the fireball 
(number of gluons per baryon, number of light quarks and antiquarks per baryon, number of 
anti-strange quarks per baryon, the pressure in the fireball, baryon density and the entropy 
per baryon) in the middle section of the table |I]. 
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After the collision has ended, for times t > 1 fm/c, but probably also < 3-5 fm/c, we relax 
the strange quarks to their equilibrium abundance and the temperature drops from Tch to 
the value Tq. The bottom portion of the table corresponds to this full chemical equilibrium 
(with exception of the S-W case for which we assume that strange quarks have reached 
80% of phase space occupancy as suggested by the experimental results 0, |38|). During 
the formation of the strangeness flavor there is already evolution of the fireball outside of 
the collision region and we allow for this by keeping Aq = Const. This effectively freezes 
the entropy content of gluons and light quarks, allowing for significant drop in pressure and 
some cooling due to conversion of energy into strangeness. 

In order to have some understanding of the possible thermal conditions prevailing in the 
early stages of the collision process, when the thermal equilibrium is reached, we have in the 
top section of the table |I] selected some reasonable off-chemical equilibrium conditions — we 
consider 20% occupancy for gluons and light quarks, and 3% for strange quarks and solve 
the same equations as for Tch, and the associated Aq . There is no change in the pressure, as 
the dynamical compression is present at this stage of the fireball evolution. But we see here 
in particular that the temperature Tth is considerably higher, since the number of quarks 
and gluons present is considerably lower. 

The columns of table |l| correspond to the cases of specific experimental interest, see top 
legend, in turn: Au-Au collisions at AGS, possible future Pb-Pb collisions at SPS with 40A 
GeV, S-Pb at 200A GeV, and for the Pb-Pb collisions at 158A GeV we considered two 
possible values of stopping, see Eq. (|2l|) : r] = 0.75 and r] = 1 . 

As discussed in section |[ the temperature values shown in the bottom portion of the table 
are similar to the inverse slopes observed in particle spectra and shown in Fig. ^. Remarkably, 
the values of temperature Tq found for the case of E/B = 8.6 GeV at ?7 = 0.5 ± 0.1 is just 
233 MeV, which corresponds nearly exactly to the reported inverse slopes of the WA85 
results [0, and Aq = 1.49 also agrees exactly with the results of our analysis 0, also shown 
in Fig.H. Even though there are a number of tacit and explicit parameters (in particular 
T] = 0.5, as = 0.6) which enter this result the degree of the agreement is stunning and 
encourages us to explore in systematic fashion the variation of the key parameters with the 
energy content of the fireball. 

In Fig.|^ we show, as function of the specific energy content E/B, the expected behavior 
of temperature Tq, the light quark fugacity Aq and entropy per baryon S/B at the time of full 
chemical equilibration in the QGP fireball. The range of the possible values as function of t] 
is indicated by showing results, for 77 = 1 (solid line), 0.5 (dot-dashed line) and 0.25 (dashed 
line). The experimental bars on the right hand side of the Fig. ^ show for high (8.8 GeV) 
energy the result of analysis of the WA85 data ||2^. The experimental bars on the left 
hand side of the Fig. ^ (2.6 GeV) are taken from our analysis of the BNL-AGS data |jT^, but 
note that in this case we had found As = 1.7 and not As = 1 as would be needed for the QGP 
interpretation at this low energy. For the BNL-AGS range of energies E/B = 2.3-2.6 GeV, 
we expect rj ~ 0.9-1 and this is indeed in good agreement with the QGP-based evaluation 
g of the experimental resuhs (T = 180 ± 30 MeV, Aq = 4.8 ± 0.4 and S/B = 13 ± 1). We 



can explain in this analysis the case of S-S collisions which have a lot of flow we 
take at E/B = 9.6 GeV a stopping fraction t] ~ 0.3 ± 0.1 which yields Tq ~ 196 ± 13 MeV, 
in agreement with the experimental results for the inverse slope [^]. The value of Aq which 
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traces the baryon density cannot be so easily estimated in this case, we refer to a study of 
the possible rapidity dependence of Aq made recently p^j . 



0.3 




23456789 10 11 



E/B [GeV] 

Figure 4: Temperature Tq, light quark fugacity Aq and entropy per baryon S/B at the time of full 
chemical equilibration as function of the QGP-fireball energy content E/B. Results for momentum 
stopping rj = 1 (solid line), 0.5 (dot-dashed line) and 0.25 (dashed line) are shown. See text for 
comparison with analysis results. 

Among the key features in the Fig.^, we note that, in qualitative terms, the drop in 
temperature with decreasing energy and stopping is intuitively as expected, and the value 
of Aq is relatively insensitive to the stopping power, and also varies little when the energy 
changes by ±15%. This implies that even when different trigger conditions lead to different 
stopping fractions rji, the resulting value of Aq which is determining the strange particle 
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(baryon/antibaryon) ratios, is rather independent of different trigger conditions. Our analy- 
sis shows that Aq decreases while E/B increases. This behavior can be argued for by noting 
that baryon density is higher in the QGP at lower energies. However, note that this intuitive 
insight was really arising from our belief that the stopping of baryon number decreases as 
energy of the collision increases, while the result here found occurs irrespective of the change 
in baryon stopping. This has considerable impact on the behavior of strange particle ratios 
as function of collision energy (see section |^). Another important result is the rapid rise of 
specific entropy with the energy content: while at the BNL-AGS energies we find similar en- 
tropy contents in the confined and deconfined phases of hadronic matter, at CERN energies 
we encounter twice as much entropy in the deconfined phase, which leads to a noticeable 



excess in particle abundances |]TT|. We turn to discuss this important issue now. 



2.7 Entropy and Particle Multiplicity 

An important observable of the thermal fireball which allows to cross check the validity of 
our approach is its entropy content. We recall the first law of thermodynamics in the form: 

E = TS - PV + i^bB , (22) 

which with the relativistic EoS (negligible masses, no dimensioned scale in the interactions): 

P = 1(1-Ab), (23) 

leads to the simple relation: 



3 W 



S AlE ^ ^ ABV 

5 = 3T5-^"^^-35T- ^2') 

Note that this equation can only be applied to systems in chemical equilibrium. The two 
last terms in Eq. (|^) are comparatively small compared to the first term. We see that the 
entropy content can be relatively accurately estimated. For E/B = 8.8 GeV and T = 235 
MeV we find S/B = 50 in agreement with other studies 0, [ll| and the results shown in 
the table above. More complete numerical calculations yield the specific entropy at fixed 
specific energy as function of one statistical parameter, or as shown in Fig.^, as function of 
the momentum stopping. In Fig. ^ we show how the entropy per baryon varies as function 
of stopping for different energies per baryon E/B = 2.6 to 10.6 GeV. The experimental bar 



in lower right corresponds to our analysis of the BNL-AGS results |]14|, and is in agreement 
with stopping being greater than 80% if QGP were formed at these energies. The cross in 
the middle of the figure is indicating the region in which we believe the S-W/Pb system to 
be (stopping 0.5 ± 0.1, specific entropy as discussed in and just above, see Eq. (0). We 
note that the solid line, corresponding to the energy content of S-W/Pb collisions, passes 
right through this region. 

This fact that entropy per baryon is primarily dependent on the specific energy content 
in the fireball offers us the opportunity to explore further the hypothesis which had lead us 
to determine the relation between the kinematic energy and the fireball energy, see section 
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E/B = 2.e 



0.2 0.4 0.6 0.8 1 

Figure 5: Specific entropy S/B as function of momentum stopping r/p at given energy content of the 
QGP-fireball E/B = 2.6, 4.6, 6.6, 8.6, and 10.6 GeV. See text for discussion of the experimental 
points. 



Should the energy stopping be greater than baryon stopping even by 15%, then the 
specific energy content would be much more greater, and in Fig. ^ we see that this would be 
implying a greater momentum stopping rjp —>■ 1. Since both energy and momentum stopping 
should increase together this is a possible scenario, and hence we must always remember that 
the 160-200A GeV collisions as possibly leading to a greater specific energy. Conversely, if 
we were to assume that energy stopping is smaller by 15% than baryon stopping, this would 
imply a reduction of the specific energy and as seen in Fig. || would require rjp 0.2, which is 
clearly not seen in the experiments. We conclude that either all stopping fractions at CERN- 
SPS energies are similar or that there is a somewhat smaller baryon stopping. The qualitative 
results about transverse energy production |^ suggest that we should not consider ?7e > 0.5 



for 200A GeV S-W/Pb collisions. 

The experimentally measurable quantity is not the entropy, but the final state particle 
multiplicity [|ri|, |l^. A useful observable is the quantity Dq: 



( d.N+ _ dN- \ 
\ dy dy ) 



\ dy dy 



Note that in the numerator of the charge of particle pairs produced cancels out and it 
is effectively a measure of the baryon number; thus in the product with S/B the baryon 
content cancels, and the result is roughly the entropy content of the final state per final 
state pion number (denominator of -Dq). Specifically, Dq ■ {S/B) ~ 3 [0] for a wide range 
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E/B 

Figure 6: Dq versus E/B. Thick line DqS/B = 4, thin line DqS/B = 3. 77 = 1 (solid line), 
T] = 0.5 (dot-dashed line) and r] = 0.25 (dashed line). 



of T, Aq and at As = 1, and taking the population of heavy hadronic resonances emitted 
from the fireball as expected from a hadronic gas system in chemical equilibrium. However, 
the hypothesis of the equilibrated hadronic g ClS dbS db transient phase is highly inconsistent 
with the strange particle flow, which suggests no re-equilibration and direct disintegration 
of the fireball into final state hadrons. If this were the case, the production of heavy meson 
resonances would be suppressed against the thermal equilibrium expectations. Consequently, 
the dilution of the number of charged hadrons by decays of heavy neutrals would not occur 
and we should observe Dq ■ (S/B) ~ 4.4, where the superscript '0' indicates that we did 
not account for the particle decay in the quantity Dq . In Fig. ^ we show as thick lines, 
using Dq ■ (S/B) ~ 4, how the observable Dq depends on collision energy. Thin lines are 
for Dq ■ {S/B) ~ 3 — the momentum stopping varies between 1 (solid lines), 0.5 (dot 
dashed) and 0.25 (dashed). Experimentally, in S-Pb collisions the EMU05 collaboration 
finds Dq = 0.085 ± 0.01, shown at 8.8 GeV in Fig.|. As the CERN-SPS energy is lowered 
to 4.3 GeV in CM frame (40A GeV for the projectile) which is presumably accompanied by 
an increase in stopping to = 1, the value of Dq rises and reaches Dq = 0.16, nearly twice 
the S-W/Pb results. We thus see that the widely different specific entropy values shown 
at the bottom of table |I| become measurable by means of the observable Dq and conclude 
that Dq is an extraordinarily valuable experimental variable which differentiates structures 
arising in the collision and from which in particular quantitative information about specific 
entropy can be gained. 
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We note in passing that a computation of the confined hadronic gas entropy p|, |TT]] at the 
same statistical conditions of the fireball yields only 50% of the QGP-fireball specific entropy 
value, and thus suggests in view of the observed particle multiplicities that the fireball could 
not have been just an assembly of confined hadrons. 



3 Strange Particles from QGP 
3.1 General Remarks 

The observable 'strangeness' is more than just one quantity which is enhanced by a factor two 
or two and a half when one compares usual p-p and p-A with A-A reactions. The interesting 
aspects of this observable are that certain strange particles are much more enhanced than 
others, and moreover, the production mechanisms being very different from the usual ones, 
the behavior of the yields (cross sections) with energy can be expected to be of different 
nature. A thermal fireball source can emit hadronic particles during its entire evolution. The 
observable strange particles are: K^*", Ks, (p, A, A, E~ , E~ , Q and Q . We should include here 
the closely related antiprotons p, which leads us to consider 11 types of particles. Leaving out 
an overall normalization factor associated with the reaction volume, and recalling that there 
is a constraint between the abundance of kaons (K^+K" ^ 2K5) we have 9 independent 
normalization parameters describing the yields. These can be redundantly measured with the 
help of the 36 = 9-8/2 independent particle yield ratios. Aside of the yield normalization 
parameters, there are in principle 11 different spectral shapes which we above implicitly 
presumed to be related to each other: The experimental fact that once effects related to 
particle decays and (transverse) flow are accounted for, the m± spectra are characterized by 
a common inverse slope parameter, cannot be taken lightly and suggest strongly that the 
source of all strange particles is indeed thermalized with a common temperature and that a 
common mechanism governs production of different strange particles. 



3.2 Final State Particle Counting 

It is easy to propose mechanisms that would largely erase memory of the transient QGP 



phase. We will not discuss this class of re- equilibrating hadronisation pictures [|2| of strange 
particles which seems not to be present at least at CERN-SPS energies 0. Instead, we shall 
focus our attention on the alternative that the particles emerge directly from the deconflned 
phase, be it because of their radiative emission be it because a general rapid and explosive 
disintegration of the hot flreball occurs in the flnal stage of its evolution. 

The abundance of particles emerging in explosive disintegration or radiated is, according 
to Eq. (^, determined by the normalization constant: 

Nj = VY[ni, ni = giXau (26) 



where it is assumed that the flnal state particle of type j contains the quark valence com- 
ponents of type i and these are counted using their statistical degeneracy gi, fugacity 
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Xi = exp{fii/T) and the chemical equihbration factor 7j. is the emission source vol- 
ume. Fragmentation of gluons could contribute to the abundance of the valence quarks and 
has been considered previously |p. Because it enhances the number of all quarks and the 
effect is weighted in a similar way for all flavors, and further, since in the ratio of parti- 
cle abundances a partial cancelation of this effect occurs, this effect is apparently of lesser 
importance. 

Once the factor 7j accounts for the deviation from the full phase space occupancy of 
the species i, the chemical potentials for particles and antiparticles are opposite to each 
other and the particle and antiparticle abundances are related, see Eq. (|]). As indicated 
in Eq. (p6D , the fugacity of each final state hadronic species is the product of the valence 
quark fugacities. The abundances of the final state strange particles is gauged by considering 
the Laplace transform of the phase space distribution of strange particles, which leads to a 
partition function Zs like expression (|27|) . The weight of individual components is controlled 
by the non-equilibrium coefficients 7s, relative meson and baryon abundance parameters C|, 
and by the fugacities Aq, Ag: 

InZ, = ^ {(AsA-i + K'X^HCl^Fj, + (A.Aj + \:'\-'HC^^Fy 

+ K'X^'HC^^F^ + {\! + K'h^C!,Fn} , (27) 

where the kaon, hyperon, cascade and omega degrees of freedom are included. Here T is 
the freeze-out temperature. The phase space factors Fi of the strange particles are (with Qi 
describing the statistical degeneracy): 

F^ = Y.9^,Wim,JT) . (28) 
j 

In the resonance sums J2j known strange hadrons should be counted. The function 
W{x) arises from the phase-space integral of the different particle distributions /(p). For 
the Boltzmann particle phase space, when the integral includes the entire momentum range, 
the well known result is found 



W{x) = (47r)-i J d\p/T)f{p) = x'K^ix) , (29) 

where x = m/T and K2{x) is the modified Bessel function. 

Because the emission volume of the particles is not known, only the ratio 

Rh = C^m/Cb (30) 

appears in observables of interest. Note that, only because we are allowing for the possibility 
of the QGP formation, we must allow for the effect of non-equilibrium hadronisation with 

M 7^ 1 — there is no reason whatsoever to expect that the rapid disintegration of 
the deconfined state will lead to particle abundances that are associated with chemical 
equilibrium of the final state. We note that relative abundances of meson and baryons 
emerging from hadronising QGP are difficult to equilibrate, because processes which convert 
meson into baryon-antibaryon pairs are relatively slow. 

Determination of the freeze-out temperature and the non-equilibrium hadronisation pa- 
rameter Rq requires that we observe the kaon to (strange) baryon ratios. We will discuss 



this point below in section 3.4 . 
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3.3 Strangeness Conservation 

There is a strong constraint between the two fugacities Aq, and As arising from the require- 
ment of strangeness conservation which was discussed at length recently 0. The manner 
in which this relation works, depends critically on the nature of the particle source and we 
shall now explore the meson-baryon off-equilibrium effects. These non-trivial relations be- 
tween the parameters characterizing the final state are in general difficult to satisfy and the 
resulting particle distributions are constrained in a way which differs considerably between 
different reaction scenarios which we have considered in detail: the rapidly disintegrating 
QGP or the equilibrated HG phase. These two alternatives differ in particular by the value 
of the strange quark chemical potential fis'- 

1. In a strangeness neutral QGP fireball Hs is always exactly zero, independent of the 
prevailing temperature and baryon density, since both s and s quarks have the same 
phase-space size. 

2. In any state consisting of locally confined hadronic clusters, /ig is generally different 
from zero at finite baryon density, in order to correct the asymmetry introduced in the 
phase-space size by a finite baryon content. 

At non-zero baryon density, that is for /ib = 3/iq 7^ 0, there is just one (or perhaps at most a 
few) special value /^b(^) which {s) = (s) at fif^ = 0, which condition mimics the QGP. 
For the case of a conventional HG (Hagedorn type) we have studied these values carefully 



10 1 . For the final state described by Eq. (UW, this condition of strangeness conservation 



takes the simple analytical form [|1T 



/x:; = Tcosh-^(i?^^-7s-^), for/z^ = 0. (31) 



There is at most one non trivial real solution for monotonous arguments of cosh~^, and only 
when this argument is greater than unity. 

Clearly, the observation of Ag = 1 (/ig = 0) is, in view of the accidental nature of this 
value in the confined phase, a rather strong indication for the direct formation of final state 
hadrons from a deconfined phase, in which this is the natural value. We note that a further 
refinement [^] of the original analysis of the S-S system at 200A GeV, which allows 
for a rapidity dependence of Aq due to flow further underpins the finding Ag = 0. We can 
thus safely conclude that strange particles produced in 200A GeV Sulphur interactions with 
diverse targets lead to a particle source which displays a symmetry in phase space size of 
strange and anti-strange particles. A natural explanation is that such a source is deconfined. 
It will be very interesting to see, if this behavior is confirmed for the Pb-Pb system and 
different collision energies. 



3.4 Freeze-out Conditions 

This is not the place to develop a complete hadronisation model of the QGP fireball — 
instead we are trying to circumvent the need for a detailed model of hadronisation by intro- 
ducing a global single parameter Rq which will allow to determine the relative ratio of meson 
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and baryon particle yields. There is a simple way to fix Rq-. the strangeness conservation 
condition Eq. ( pl| ) must apply to the particle emitted from the fireball, with chemical condi- 
tions determined by the source: thus for each value of /iq and T = Tf there is a corresponding 
value of Rq for which the strangeness remains balanced. Rq = 1 would almost always lead 
to asymmetric emission of the strange particles and thus in all cases to the formation of a 
final s or s strangelett nuggets, except in the unusual condition that the QGP disintegration 
occurs when the plasma and HG phase-spaces for s or s particle-carriers are the same. In 
any event, it is extremely instructive to study how the values of the so defined Rq depend 
on the properties of the fireball. 

We use here the experimentally motivated 7s = 0.7, though the deviation from unity is 
of little numerical importance in present argument. Since the (multi-) strange (anti)baryon 
particle ratios have led to Aq = 1.48 and As = 1.03 we use the QGP freeze-out value As = 1 
and set Aq to three values in Fig.|^. The solid line is for Aq = 1.5, an appropriate choice for 
the case of S-W collisions at 200A GeV, (when Ag = 1). The short dashed curve is for the 
choice Aq = 1.6 which is as large as Aq will get for the case of Pb-Pb 160A GeV collisions; 
we see that we are finding rather narrow range of Rq for each given freeze-out temperature. 
The long dashed curve is for Aq = 2.5, the value which our model calculations suggest for 
the 40A GeV collisions (see table |I|). 

The value Rq = 1 is found for T ~ 200 MeV at Aq ^ 1.5-1.6. For lower temperatures 
we have Rq < 1, which implies that the size of the strange baryon phase space needs to be 
increased compared to thermal equilibrium (or the strange meson phase space needs to be 
reduced) in order to allow a balance between strange and anti-strange quarks in emission 
from the QGP fireball. 

As noted in passing above, a description of hadronisation using Rq = 1 leads to a 
s, s asymmetry in the developing fireball. The question is which of the two evolution 
options: the off-equilibrium particle emission, or strangeness distillation, is consistent with 
the experimental data. This can be decided considering a physical observable which is 
primarily sensitive to the parameter Rq, and to a lesser degree to the other thermal model 
parameters. We choose to study the kaon to hyperon abundance ratio, both for all phase 
space and at fixed m±. In the latter case we have: 



R I - 



A + SO 



As/ Ad + Ad/ As ^22^ 
8 AgAuAd 



We made a preliminary study of this relation in |TD[, which is valid when resonance decay 



contributions cancel. To account here in necessary detail for the resonance decay influence, 
we have incorporated the decay pattern of all listed resonances numerically. The thick lines 
in Fig.Hgive -RkIttix fixed m± as function of Rq, see Eq. (^2]). We include the descendants 
of strong and weak decays in order to facilitate comparison below with experimental data. 
We assumed that the distribution of parent particles for kaons and hyperons is according 
to the thermal equilibrium condition evaluated at temperature as given in the Fig.|^. The 
values of Aq = 1.5, 1.6, 2.5 and line conventions are the same as used in that figure. The 
thin lines in Fig. ^ give results covering the full phase space. Since the kaon mass is much 
smaller than hyperon mass, this ratio -RKltot is considerably greater than -RKlmx; hence 
we use logarithmic scale to display both i?K- We see in particular that in order to have 
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Figure 7: Strangeness neutrality line: Rfj versus freeze-out temperature for Aq = 1.5 (solid line), 
Aq = 1.6 (long dashed line) and Aq = 2.5 (short dashed line). 



HG-equilibrated yield of kaons and hyperons for fixed m±, i.e. when Rq = 1, we would have 
to find experimentally -Rxlmx — 0-3 in the S-W/Pb and/or Pb-Pb collisions. 

There is no officially reported value for the i?K ratio. However, collaboration WA85 has 
presented [^] a figure which shows for largely overlapping 1.1 < m± < 2.6 GeV the yields 
of A, A and Kg for the central rapidity region 2.5 < y < 3. No cascading corrections was 
applied to these results. Using graphic methods we obtained -Rxlmx = O-H i 0.02. This 
implies a far off-HG-equilibrium result Rq = 0.4 as shown in Fig.^, and thus as seen in 
Fig. 0a freeze-out temperature Tf ~ 150 MeV. The possibility that Rq ~ 1 (-Rxlmx — 0.3) is 
experimentally completely excluded. The factor Rq ^ 1 confirms the expectation that the 
reactions which change the number of baryons (baryon-antibaryon formation by mesons) are 
relatively slow in the confined phase , or that all strange particles originate directly from 
the QGP fireball. In any case we conclude that the observed value Rq = 0.4 allows to fix 
the freeze-out conditions, here in particular Tf = 150 MeV. 

The final issue is how, from the value Rq = 0.4, we can infer the values of the abundance 
constants Cf^ and which (see Eqs. (p7|,|30D) express the relative strange meson and baryon 
production abundance to the thermal equilibrium values. If we argue that the strange meson 
abundance, akin to total meson abundance is enhanced by factor two (i.e. = 2) as we 
found studying the entropy enhancement, then the conclusion would be that the strange 
baryons are enhanced (against their tiny HG equilibrium abundance at Tf = 150 MeV) by 
factor = 5. It is clear from this observation that strange baryons and antibaryons are 
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Figure 8: Thick lines: RK\mj_, thin hnes: i?K|tot as function of Rq for Aq = 1.5 (sohd hne), Aq = 1.6 
(long dashed line) and Aq = 2.5 (short dashed line). 

thus the key to the study of the deconfined state. We now turn to study these particle yields. 
3.5 Strange Baryons and Antibaryons 

The ratios of strange antibaryons to strange baryons of same particle type: Ra = A/A, 
i?s = S/S and i?n = ^/^, are in our approach simple functions of the quark fugacities and 
have been discussed at length recently [§, |^. In view of our current work we can predict 
the behavior of these ratios as function of energy. Using the results for Aq shown in Fig.^, 
we show in Fig. |] these three ratios. Since we assume As = 1 , we have = = 1, but 
since some re-equilibration is to be expected towards the HG behavior As > 1, we expect 
As = 1 + e, with e small, and thus for this ratio Rq = 1 — 6e < 1. A further non negligible 
correction which has been discussed at length in Ref. [0] is due to the isospin asymmetry: in 
the heavy Pb-Pb collisions it will be necessary to account for d-u asymmetry which is as 
large as 15%, and which favors the abundance of particles with rf-quark content over those 
with M-quark content. This impacts here in particular the ratio R^, since there are no light 
quarks contributing to Rq, and the ratio Ra is u-d symmetric. 

We now study the ratios between antibaryons with different strange quark content. These 
are dependent on the degree of the strangeness saturation reached, and we shall take in this 
study 7s = 1 assuming relatively large, long-lived system created in the collisions of largest 
available nuclei. In the Figs.p!0|-[12| we show three ratios and for each ratio three results: 
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Figure 9: Antibaryon to baryon abundance ratios as function of energy per baryon E/B in a 
QGP-fireball: i?N = p/p (solid line), R\ = A/A (long-dashed line), i?s = H/H (short-dashed line) 
and Rq = 0/17 (dotted line) 



solid lines depicts the result for the full phase space coverage, short dashed line for particles 
with p± > 1 GeV and long dashed line for particles with m± > 1.7 GeV. In Fig. ^ we show 
the ratio A/p, in Fig.[TT| the ratio S~/A and in Fig. [12| the ratio Q/E^. Because Aq rises with 
decreasing E/B, see Fig.^, we encounter overall the remarkable behavior that these three 
ratios increase as the collision energy is reduced. 

The behavior shown in Figs. p!0|-[T^ may be of considerable importance, if our expectations 
are confirmed that in reaction models in which QGP is not assumed and the particles are 
made in a sequence of microscopic collisions these ratios do increase with the collision energy, 
reflecting in this behavior the behavior of the reaction cross section observed in p-p reactions. 
If this behavior is found at low energies and it matches at some energy at which a jump to 
the here presented yields, then the QGP behavior obtained here is truly the smoking gun 
type evidence for the formation of the deconfined phase. Furthermore, it would be a rather 
easy task to determine the transition energy to QGP by merely seeking where these ratios 
peak as function of E/B. 

As a final step in this discussion we present now the analysis of the available and very re- 
cent WA85 fi/S~ production ratio and the K/p ratio of the NA35 collaboration obtained 
for the S-Au system at 200A GeV [^]. Fig. shows a comparison of our ab initio calculation 
and the pertinent experimental results. We use the same cuts on the range of p_[_ as in the 
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Figure 10: Strange antibaryon ratio A/p, as function of E/B in a QGP-fireball; solid lines are for 
full phase space coverage, short dashed line for particles with p± > 1 GeV and long dashed line for 
particles with mj_ > 1.7 GeV. 
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Figure 12: Strange antibaryon ratio il/H^ with the same conventions as in Fig.|lC 
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Figure 13: Strange antibaryon ratios for S-W/Pb collisions as function of E/B in a QGP-fireball: 
A/p (full phase space), E= /A for p_L > 1.2 GeV and (H + n)/(E^ + H-) for px > 1.6 GeV; 
experimental results shown are from experiments NA35, WA85. 
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experiment: the experimental points show the results A/p ~ 0.8 ±0.25 (NA35) for full phase 
space, S^/A = 0.21 ± 0.02 (WA85) for > 1.2 GeV; and (fi + n)/(S- + H^) = 0.8 ± 0.4 
(WA85) for > 1.6 GeV. The chosen value of 7s = 0.70 and rj^ = 0.5 brings about good 
agreement of our model with the precise value of S~/A. Fig. ^ shows also the impact of the 
change of the collision energy on these results, using 50% stopping, rather than rj = 1 used 
in Figs. |10H1|. 

Considering that we have computed here everything in an ab initio dynamical model 
(which as discussed above has some tacit and explicit parameters such as the QCD coupling 
as = 0.6 etc., chosen to be in agreement with the earlier experimental results) it is remarkable 
that such a good agreement with the two very recent results could be attained. We can 
conclude that the fact that the two ratio A/p (NA35) and {Q + + E^) (WA85) are 

satisfactorily explained, provides a very nice confirmation of the consistency of the thermal 
fireball model. 

An interesting question which arises quite often is how the particle and in particular 
antibaryon yields vary with energy. Eq. (^) allows to determine the absolute particle yields 
as function of fireball energy. Considerable uncertainty is arising from the off-equilibrium 
nature of the hadronisation process, which in particular makes it hard to estimate how the 
different heavy particle resonances are populated, and also, how the abundance factors Cg 
vary as function of energy. Some of this uncertainties are eliminated when we normalize 
the yields at an energy, which we take here to be the value E/B = 2.6 GeV which is 
applicable to the BNL-AGS experiments. In Fig. |14| the so normalized yields of antibaryons 
taking the freeze-out temperature T = 150 MeV are shown (we also assume 7s = 1, //p = 1 
and absence of any re-equilibration after particle emission/production). These yields are 
rising in qualitatively similar systematic fashion with energy, as would be expected from 
the microscopic considerations, but the rise of more strange antibaryons is less pronounced, 
unlike what we would naively have expected. The quantitative point to note is that at BNL- 
AGS {E/B = 2.6 GeV) the yield from a disintegrating QGP-fireball is a factor 100-400 
smaller compared to yields aX E/B =9 GeV. Since the particle density dN/dy is not that 
much smaller at the lower energies (recall that the specific entropy, see table |I| , drops only 
by factor 3.5, implying a reduction in specific multiplicity by a factor 5), it is considerably 
more difficult at the lower energies to search for antibaryons than it is at higher energies. 
We should remember that the results presented in Fig. |1^ are obtained assuming formation 
of the QGP-fireball and same freeze-out and hadronisation conditions for all energies shown. 



4 Summary and Conclusions 

The key issue of interest to us is the identification of the deconfined quark-gluon plasma 
state in relativistic nuclear collisions. We have shown in recent work that many features 
of strange particle production results obtained at 200A GeV, are consistent with the QGP 
hypothesis; and we have argued that the conventional reaction picture involving a fireball 
made of confined hadrons is for a number of reasons incompatible with these experimental 
results. We felt that in order to ascertain the possibility that indeed the QGP phase is 
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Figure 14: Relative antibaryon yields as function of E/B in a QGP-fireball. p (solid line), A 
(long-dashed line) (short-dashed line) and O, (dotted line), all normalized to their respective 
yields at E/B = 2.6 GeV . 

already formed at 200A GeV a more systematic exploration as function of collision energy 
of these observables would be needed — conclusions drawn from one set of experimental 
conditions suffer from the possibility that some coincidental and unknown features in the 
reaction mechanisms could lead to the observed properties. It is highly unlikely that this 
would remain the case, should a key feature of the experiment be varied. 

Variation of the number of participating nucleons has in the QGP reaction picture the 
considerable disadvantage that we reduce the size of the deconfined region, which even in 
the case of Pb-Pb reactions, is not very large. This very likely will reduce the hfespan 
of the deconfined state, and the key reaction measure, the stopping 'friction' is reduced in 
such a case as well. All these effects make an interpretation of the data with varying the 
number of participants, in terms of the same reaction model, exceedingly difficult, if not 
impossible. Against this drawback, the study of the particle production with change in 
the reaction energy (measurement of excitation functions) appears to be the best approach, 
which preserves the maximum attainable size and lifespan of the hypothetical deconfined 
region, while minimizing the change in the stopping, in particular if one is operating in 
the energy domain in which for the largest available nuclei the stopping is nearly complete. 
Furthermore, at some low energy, we must reach the conditions of the normal hadron cascade 
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interaction, and thus be able to describe the results reliably in terms of established cross 
sections. 

Thus our primary objective in this work was to establish the systematic behavior of 
the antibaryon yields as the energy and stopping of the colliding large nuclei vary, and 
to determine how the freeze-out conditions of these particles can be determined. This is 
necessary since the experimental results suggest that while the thermal (kinetic) equilibrium 
is established, the chemical (particle abundance) equilibrium in the processes governing final 
state particle freeze-out is largely not achieved. We thus developed a model which relates 
the energy content of the fireball to the collision energy and the stopping friction. In order 
to relate the energy to the statistical properties of the fireball which in our approach controls 
the particle yields we had to use EoS of the deconfined matter, and in particular to account 
for the impact of QCD-interactions on the properties of the QGP, see section |2]^. We have 
then shown how the study of kaon to hyperon ratios allows to determine the meson-baryon 
chemical non-equilibrium parameter. We found that while the entropy production suggests 
a meson yield which is about twice that of thermal freeze-out yield, the baryon abundance 
is five times greater than the small thermal freeze-out abundance. 

We believe that the thermal equilibrium approach should not be advanced in the context 
of A-A collisions as being just an economical, but otherwise approximate and highly limited 
method. Actually, one can argue that this is a more appropriate approach compared to 
microscopic dynamical (cascade) models, pending a better understanding of the elementary 
soft particle and entropy producing hadronic processes. In the A-A reactions there is no 
reason to expect that 'cascades' of p-p type interactions, which are, as matter of principle, 
even less understood, lead to adequate understanding of soft particles in nuclear collisions. 

We have used a simple model which allows to determine, in a systematic fashion, the ther- 
mal conditions reached in high density deconfined matter generated in heavy ion collisions. 
It is based on the observation that during the collision the compression of the quark-gluon 
matter can proceed until the internal pressure succeeds in stopping further impact com- 
pression. In this picture we worked far from chemical quark-gluon equilibrium, but always 
assumed that the thermal equilibrium is reached rapidly, even on the scale of the nuclear 
collision time, which is believed to be ~ 1 fm/c given the geometric size of the nuclei and 
the Lorenz contraction factors. 

Motivated by the absence of chemical particle abundance equilibrium, natural for a ther- 
mal HG in the final state, we employ a picture of particle production which involves rapid 
disintegration of the QGP-fireball . Central to the particle abundances are then the chem- 
ical properties of the QGP-fireball and we have discussed these comprehensively as function 
of collision energy and stopping in section |2.6| . We have shown that the thermal conditions 
we find at the end of strangeness chemical equilibration in the fireball, see bottom of table 
0, are in good agreement with our expectations derived from particle yields seen in S-Pb/W 
collisions. To carry through this program we needed to make a reasonable choice of remaining 
parameters: in particular stopping rj = 50% — about equal for baryon number, energy and 
pressure; perturbative EoS of particular form we discussed at length in section Given 
these assumptions and parameters, we were able to study, the current strange particle data 
at 200A GeV and have reached a very satisfactory agreement with experiment as is shown 
in Fig. [To. 
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With this encouragement, we have computed in a systematic fashion the behavior of 
strange particle (kaon, baryon and antibaryon) yields assuming conditions likely to occur in 
Pb-Pb interactions (e.g. full stopping and 7s = 1). It is most interesting that these results 
show patterns of behavior which could indeed be unique for the QGP type of fireballs — in 
particular, the relative yields of strange antibaryons (see Figs. ^ - [12|) lead to greater ratios 
occurring at smaller collision energies, down to very low energies, where in models involving 
cascading type reactions such ratios would be so small that their measurement would be 
difficult. We are persuaded that this pattern of behavior could not occur for normal confined 
matter, where the rise in cross sections with energy dominates particle yields whenever these 
are not arising from collective phenomena such as is a deconfined QGP phase. 

A much discussed question p|, has been whether current strange particle data are 

consistent with the fully equilibrated hadronic gas (HG) picture of hot hadronic matter. The 
motivation for such a hypothesis emanates from the observed rapid thermalization: given 
that we do not understand the mechanisms of thermalization, one could equally argue for the 
'a priori' presence of fully equilibrated (thermal and chemical) HG phase. We do not believe 
that the chemical and thermal equilibration are due to same physical processes, and thus 
while allowing thermal equilibrium we maintain the possibility in our analysis that chemical 
equilibrium develops slowly and not completely during the reaction. The results presented 
here support this point of view. 

We stress that our description of particle production is based on collective mechanisms 
(QGP-fireball) and is thus intrinsically different from microscopic approaches, in particular 
when these are based on a hadronic cascade picture. Such models generally exploit specific 
data and/or extrapolations and assumptions about individual hadronic reactions and their 
cross sections. If the true underlying thermalization processes are different (as is likely) from 
those used in current microscopic approaches, and are indeed much faster than cascades 
suggest, the whole representation of the collision evolution in these microscopic approaches 
needs to be reviewed — in this perspective our thermal approach, already advocated for 
the p-p reactions maybe seen as an experimentally well motivated hypothesis. We 
note that no alternative model to the here developed rapidly hadronising QGP has been 
proposed which could generate both strangeness abundance and multi-strange antibaryon 
enhancement. For example the description in terms of the dual parton model DPM, | ^ , 
which introduces a number of parameters to enhance strangeness and strange antibaryons, 
arrives at a considerably smaller relative abundance . 

Some of our results, though presented in great numerical detail, could suffer from consid- 
erable systematic errors, in particular when overall absolute yield normalization is discussed. 
For example, to determine the freeze-out conditions we need to determine the kaon yield af- 
ter all heavy resonances have decayed. We employed all known tabulated resonances (and 
hence probably not all), and we adopted the statistical spin/isospin factors not always es- 
tablished. Furthermore, there is the influence of the large resonance width which allows that 
these states are formed at energy below the resonance mass. There is also the possibility 
that some resonances alter their properties (particle width, mass) in dense matter. All these 
effects may have signiflcant impact on the relation we have established between kaon to 
hyperon ratio and thus on the freeze-out conditions, and off-equilibrium properties, which 
in turn impacts the individual particle yield. The presence of these effects and uncertainties 
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is also behind our omission from this discussion of the 0(ss)-meson, which participates and 
benefits from the general strangeness enhancement effects, but which is difficult to describe 
with the precision that is needed. 

Fortunately, these systematic uncertainties have very small impact on the strange an- 
tibaryon ratios, in particular when these are considered at fixed, high m±, and we firmly 
believe, in view of the results we have obtained here, that such data provide the best hadronic 
signatures, and diagnostic tools, of the deconfined matter. We recall the large ratios in the 
QGP-fireball reaction picture, such as S/A which we have found at relatively small energies 
— in microscopic models and near to S production threshold in p-p interaction this ratio is 
very small. This lets us expect that there will be a peak in the relative S/A yield as func- 
tion of collision energy which will provide an interesting possibility to identify the energy at 
which collective production of strange antibaryons is first encountered. At this energy we 
should also encounter for the first time the other features of the QGP phase: strangeness 
production enhancement, strange phase space saturation (7s — > 1) , entropy enhancement 
(particle multiplicity enhancement), pattern of strange antibaryon flow showing Ag = 1. 
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